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1 In a recent Letter Jy, Geiger presents calculations of the dilepton emission from the 

m : 

early stage of ultrarelativistic heavy ion collisions using the parton cascade model (PCM). 

O; 

He shows that the M± scaling found in Ref. is not observed. In this Comment, we 
^ | point out that this is largely due to a defect in the PCM. 

■ The PCM is based on the perturbative QCD and inevitably contains two infrared 

cutoffs. In particular, it has the lower limit of the timelike virtuality of partons, /io = 
lGeV. As a result, each parton has an effective mass larger than lGeV. During the nuclear 
collision, many partons acquire large timelike virtuality and it decreases gradually to /i . 
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In the PCM, once the virtuality reaches /io, it is kept unchanged. This is expected to 
happen typically by r ~ l//i = 0.2fm = r cut , where the proper time r is defined as 



r = sgn(t-to) y (t — to) 2 — z 1 [yj and in the center of mass frame, at z — 0, the maximum 
density is achieved at r = 0. In reality, the virtuality of partons continues to decrease 
below /io as time increases according to the uncertainty principle. Therefore, the partons 
in the PCM have inappropriately large virtualities after r ~ r cut . 

This unphysical virtuality drastically affects the dilepton production, especially that 
from qq annihilations with M comparable to ~ 2/io or less, where M is the invariant 
mass of the dilepton. A remarkable feature in Fig. 3 of Ref. Q, on which the formalism 
of Ref. m is based, is strong suppression of the dilepton yield from qq annihilations 



below M ~ 2.8GeV. At a glance, this might look like a non-equilibrium effect. This 
is, however, not the case. Note that schematic calculations of the dilepton emission 
from the preequilibrium stage with on-shell quarks show enhancement of the dilepton 
production from the qq process at small M [|5]. Since this is a QED process and there 
is no divergence in the process, the other cutoff in the PCM, p± cu t, does not play a role 
here. This suppression is due to phase space suppression by the large effective quark 
mass. In this PCM calculation, there is no qq contribution below M = 2GeV, whereas in 
conventional calculations qq annihilations dominate this region. From Fig. 3 of Ref . || , 
it is natural to expect that the qq contribution would be dominant also at M ^ 3GeV if 
it were not excessively suppressed by the large unphysical effective quark mass. 

According to Fig. 3 of Ref. [|lj , most of the Mj_ scaling breaking dileptons are created 
after r = 0.5fm, when the PCM has too large effective quark masses. In the PCM local 
thermalization is achieved as early as r ~ 0.3fm || and the transverse expansion does not 
establish in the early stage of the collisions. Therefore, neither the preequilibrium state 
nor transverse expansion is the reason for the M±_ scaling breaking. It is the dileptons from 
the bremsstrahlung that are responsible for the Mj_ scaling breaking. The bremsstrahlung 
contribution is dominant at small M in the PCM jp and does not generally realize the 
Mj_ scaling. For the M± scaling, the Boltzmann approximation is essential @, in which 
diagrams with outgoing (anti) quarks or gluons are neglected. It cannot be used for the 
bremsstrahlung. 

Since, as aforesaid, the other conditions for the Mj_ scaling are approximately satisfied 
as early as r ~ 0.3fm, if the virtuality of the partons is appropriately treated and the 
unphysical artificial suppression of the qq process at small M is removed, the large M± 
scaling breaking obtained in Ref. [JJ is expected to disappear. 

In summary, the large Mj_ scaling breaking reported in Ref. [|]] is due to the use of the 
inappropriately large virtualities > fio after r ~ r cut and the consequent suppression of qq 
annihilations at low M, and is not physical. It is not justified to use the present version 
of the PCM to calculate the dilepton yield and confirm the M± scaling in the region of 
M, ~ 2 - 3GeV. 
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